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Abstract 

Bound states in quantum mechanics must almost always be constructed numerically. 
One of the best known exceptions concerns the central D— dimensional (often called 
"anharmonic" ) Hamiltonian H = + a |r|^ + h \r\^ + . . . + z |r|^5+2 (^-yyj^gj^g z = 1) 
with a complete and elementary solvability at g = (central harmonic oscillator, no 
free parameters) and with an incomplete, A^— level elementary analytic solvability 
at g = 1 (so called "quasi-exact" sextic oscillator containing one free parameter). 
In the limit D ^ oo, numerical experiments revealed recently a highly unexpected 
existence of a new broad class of the g— parametric quasi-exact solutions at the next 
integers g = 2,3,4 and g = 5. Here we show how a systematic construction of the 
latter, "privileged" D ^ 1 exact bound states may be extended to much higher gs 
(meaning an enhanced flexibility of the shape of the force) at a cost of narrowing the 
set of wavefunctions (with restricted to the first few non- negative integers). At 
g = 4K + 3 we conjecture a closed formula for the A^ = 3 solution at all K. 
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1 Introduction 



In quantum mechanics, bound states of a particle confined in a central potential well 
y(|r|) in D dimensions are constructed as normalizable solutions of the ordinary 
differential Schrodinger equation 
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Usually, this equation (or, more precisely, their infinite set numbered by integer 
argument L = 0, 1, ... of the "angular momentum" i — £{L) = L + {D — 3)/2) must 
be solved by purely numerical means. 

A notable exception concerns all the models where the spatial dimension D 
proves "sufficiently" large, D ^ 1. In a way outlined, say, in reviews [1], a semi- 
analytic, perturbative construction of the solutions may be then based on the so 
called large— expansion technique. In essence, this approach combines our knowl- 
edge of the asymptotic growth of V{r) ~ r" at r ^ 1 with the presence of the strong 
repulsive core i{i-\-l)/r'^ near the origin. This implies the existence of a (presumably, 
pronounced) absolute minimum of the combined forces at a point r = R{a) where 
we have 
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whenever we neglect all the less relevant corrections. The main merit of this ap- 
proach lies in its simplicity and universality but, unfortunately, its expected rate of 
convergence (measured by the smallncss of 1/R) will, obviously, increase quite slowly 
with the growth of D and decrease quite quickly with an increase of the dominant 
power of the potential a. As long as we are going to pay attention to the class of 
potentials 
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with a fixed k (say, for the sake of simplicity, k — 0) and growing integers q, an 
efficient use of the above large— D approach does not seem too promising. 

In an alternative semi-analytic approach to the problem (1) -|- (2) one could try 
to employ a power-series method where 



/"^^ X U{r) X expVF(r) , 
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and where approximations U{r) — polynomial [r^/(^+^)] (of a sufficiently large degree 
mu) and W{r) — polynomial [r^/('^+^)] (of any degree mw) may be re-constructed 
in a more or less algebraic manner after an insertion of this ansatz in the original 
differential equation. Unfortunately, virtually all practical implementations of this 
method prove even less efficient in computations [2]. 

The third eligible non-numerical approach to the construction represents a modi- 
fication of the power-series method where also the second polynomial degree mu < oo 
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is fixed. This method has been proposed by E. Magyari [3] and is based on the eval- 
uation of the solutions (3) via an ad hoc tuning of the potential (2) itself. Basically, 
the recipe requires the absence of any errors in eq. (3). This idea just extends the 
well known polynomial solvability of eq. (1) + {2) a.t q — k — (harmonic oscillator) 
and q — k — 1 — (Coulomb field). The details of the Magyari's recipe have been 
reviewed in rcf. [4] where we emphasized that its practical merits are virtually non- 
existent because the underlying and obvious selfconsistency requirements (wc called 
them there the Magyari-Schrodinger equations) may be characterized as a coupled 
set of determinantal equations which are extremely complicated to solve in general. 

The main reason why we decided to pay attention to the Magyari's approach is 
that wc discovered, many year ago [5] . that this method encounters enormous simpli- 
fications during a limiting transition to the large dimensions D ^ 1. In this setting 
one has to follow very consequently all the analogies with the harmonic oscillators. 
Firstly, one assumes that the polynomial W{r) in the exponential represents an exact 
r ^ 1 asymptotic solution at any q and k. Both the degree rriw and all the coeffi- 
cients in W(r) are uniquely determined, in this way, by our choice of the potential 
(for example, we have W{r) = — r^/2 for harmonic oscillator). On this background, 
the Magyari's second key idea parallels the essence of the quasi-exact solvability (i.e., 
an incomplete solvability, see [6]) and extends, to all the q> I cases, the requirement 
that our choice of the polynomial U (r) of any integer degree N — mjj leads to some 
exact solution (3) for a suitably modified, adapted polynomial potential (2). 

Our present text is inspired by the series of demonstrations that the latter com- 
bined method proves unexpectedly efficient in practice. As wc mentioned, our study 
was initiated by the observation that for one of the most popular, viz., quartic 
polynomial interaction, the complicated Magyari-Schrodinger equations exhibit a re- 
markable and fairly surprizing simplification. In the next stage of our work on this 
project [7] we choose the "first nontrivial" q = 2 oscillators and complemented the 
above observations by the formulation of the method of an explicit perturbative con- 
struction of 1/D corrections. We must emphasize that the improved convergence 
of our innovated solutions (which were defined by the series in the powers of 1/-D^) 
was in a sharp contrast to the steady worsening of the performance of the above- 
mentioned large— D expansions at the larger q. Moreover, in the "minimal" sextic 
example at g = 1, we succeeded in showing that and why our innovated series in the 
powers of l/i? was absolutely convergent [8]. 

A real climax of our effort came with the papers [4] and [9] where we performed 
an explicit construction of the zero-order solutions and discovered that at the next 
few less trivial exponents g = 3, g = 4 and g = 5 we were still able to construct 
the D — > oo solutions in closed form. Unfortunately, we were just able to work 
by the brute-force methods, based on the direct solution of the D ^ oo limit of 
the Magyari-Schrodinger coupled algebraic nonlinear equations by the elimination 
method using the Groebner bases [10]. For this reason, we were never able to find 
any solution at g > 6. This was also the main motivation for the study which we are 
going to describe in what follows. 
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2 A concise formulation of the problem 



2.1 Asymptotic Magyari-Schrodinger equations 



In a way discussed thoroughly in [4], the use of our ansatzs in the hmit D ^ oo 
imphes an immediate reduction of the differential radial Schrodinger equation to its 
algebraic equivalent 



/ si 

S2 



N -1 



1 



\ 



N-2 Sn 



Sl N-2 

Sl N-1 

: Sl 



'9-1 



V 



/ Po \ 

Pi 
PN-2 

V Pn-1 J 



0. 



(4) 



Here we are going to study its solvability, referring to [4] (we shall call it paper PI 
in what follows) for all the explanations of its origin and interpretations. Calling 
this deeply nonlinear non-square-matrix algebraic problem the Magyari-Schrodinger 
equation, we may just summarize that the N quantities pj are in a one-to-one corre- 
spondence with the Taylor coefficients in the parts U{r) of wavefunctions (3). More- 
over, all q generalized eigenvalues Sa [cf. their definitions Sa = Sa{ga-2) by eq. (30) 
in PI] arc just certain re-scaled forms of the energies and/or of coupling constants in 
our potential (2). 



2.2 Upside-down symmetry 

As we already explained in paper PI, the practical use of the explicit quasi-exact (QE) 
solutions requires a purely numerical determination of their QE-compatible eigenval- 
ues ga-2 with a — 1,2, ... ,q. This means that the simplified Magyari-Schrodinger 
eq. (4) should be interpreted as an implicit definition of a QE solution on a pre- 
determined level of a finite precision. All the sufficiently small 0{1/D) variations of 
the parameters should be ignored, within this convention, as negligible and irrele- 
vant. Alternatively, they may all be treated, if necessary, by perturbation methods 
in a way exemplified in [11] at q — 1 and in [7] at q — 2. 

In such a setting, one of the main messages delivered by paper PI was that all 
the perturbation constructions of any type may remain exact and non-numerical, 
at the first few smallest q at least. This follows from the observation that all the 
exact QE parameters in PI proved expressible through integers in the zero-order limit 
D ^ oo. The necessary condition of this simplification lies in the elementary form of 
our algebraized QE Magyari-Schrodinger equation (4). Still, one has to overcome a 
few further obstacles. In particular, the brute-force origin of the results in PI made 
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it impossible to move beyond q < 5. In the other words, the key weakness of paper 
PI lies in the too rapid growth of the difficulties with the increasing q. One needs 
all the capacity of the available computers to reveal the structure of solutions at the 
first few q. This means that whenever one needs an improvement of the insight in 
the structure of the solutions, one has to make a better use of their symmetries. 

In our present continuation and completion of paper PI we arc going to exploit 
the most obvious symmetry of eq. (4) with respect to its upside-down transposi- 
tion. For this purpose we may modify our notation slightly, replacing eq. (4) by its 
reincarnation 
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a2 Oil 2 

Oil '■ ai N — 2 

N -I Oil ai N -I 

N — 2 OLi : «! 

2 Oil 01-2 

\ I oil J 

However trivial, such a change of notation implies that all the separate lines have a 
tilded partner in this set. For example, the definition po = —pi/oi (where a — ai 
is tacitly assumed non-zero) is always accompanied by its tilded counterpart po = 
— pi/ OL whenever a = ai 7^ 0, etc. 
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3 Solutions of the symmetrized recurrences 

3.1 Trivial case at iV = 2 

In our tilded notation, the ffist nontrivial version of our problem (5) with = 2 is 
particularly transparent and instructive, 
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Its upside-down or "tilding" symmetry separates its first and last line as giving a 
constraint upon the doublet of unknowns ai — si — a and oii — Sq — a m. ^ 
tilding-symmetric manner. 
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The fully tilding- symmetric way of dealing with the rest of eq. (6) consists now in 
the equivalence of its recurrent downward or upward treatment. 

After we normalize pi = p = 1, the first line defines po = P = — l/a and may 
be omitted. Step-by-step, the k—th line of eq. (6) may be multiplied by the factor 
a/a^~^. This converts the old components into the known constants and we get the 
following sequence of definitions 

Ota = 1, 01012I01 = 1, aa^/a^ = 1, . . . , (8) 

with elementary consequence: = Oi^. We may imagine that the last definition 
prescribes that ctg+i = a'^^^ = 1. This equation may be read as a boundary condition 
for our recurrences, fixing the physical value of our single free parameter a. It has 
many unphysical complex roots and just a single real one, viz., the physical root 
a = 1 at any even q. Similarly, two different real roots a = ±1 become available at 
all the odd g's. This makes the final reconstruction of all the original Q&compatiblity 
"eigenvalues" Si, . . . ,Sg trivial. 

We may add a comment. Knowing that the last line of recurrences (6) defines the 
function of a (ctg+i = a^"*"^) with a prescribed value (a*"'"^ = 1), the latter constraint 
may be interpreted as an algebraic equation which fixes the eligible values of a. Such 
a type of the boundary condition is not unique. The same role may be played by 
any other line of eq. (6), once we re-direct these recurrences and demand that 

ai+j = ftg-j (9) 

at any shift j < q — 1. At A?" = 2 all this is trivial since after being multiplied by 
a^'^^, all relations (9) degenerate to the same rule 0;^+^ = 1. 



3.2 A transition to the single variable at N = 3 

A separate treatment of the first nontrivial N — 3 version of eq. (4) is necessary in 
the degenerate case with pi = 0. We may infer that Si = S3 = = . . . = Sq = 0. 
This means that q = 2Q + 1 must be odd and that we in effect return to the N = 2 
structure. We only have to replace the old unknowns by the new ones, re-scaled to 
S2fe/2. Otherwise, the construction of the solutions remains strictly the same, giving 
the nontrivial roots S2k — 2^ where q^^^ — 1. 

In what follows, similar detailed qualification will be omitted and, with the de- 
generate solutions ignored, we shall normalize pi = 1 at A^ = 3, etc. 

From the two outer lines of the A^ = 3 version of eq. (4) we deduce that po = — l/a 
while P2 = — l/cK. The rest of equation (4) acquires the tilding- symmetric matrix 
form 
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It gets facilitated when pre-multiplied by an auxiliary row vector. This observation 
results from the step-by-step analysis of this system of equations re-written in the 
form 

a 13/a = aa — 2, I a = afi — a.^ a5/a = Of-y — {3, .... (11) 

In the first item the right-hand-side part aa — 2 = aa — 2 = ^ — 2 = ,^ — 2 is tilding 
symmetric. This means that the same tilding-invariance must hold for the left-hand- 
side expression as well. The second item is not tilding-invariant but the invariance 
is restored after we divide all this equation by a. This gives a consistent picture be- 
cause one can deduce that also in all the subsequent rows the full tilding-invariance 
is achieved when we replace a, (3, 7, ... by their renormalized and tilding-invariant 
forms aa/a'^, (3a/a, ^ya/a^, . . ., respectively. In the other words, the system (11) 
must be pre-multiplied by the row of the factors 1, l/a, a/a, a/a"^, 0? jo? , . . . ob- 
tained, in recurrent manner, by the multiplication by the quotient which depends 
on the parity, i.e., equals to l/a and to a in subsequent steps. This means that the 
even and odd items in eq. (11) have a different structure. 

This difference may be reflected by the change of the notation. Once we put 
a — S\ — a, (5 — S2 — A, ^ — s-i — h, ^ — s/^ = B , e = = c (while denoting also 
Sq — a etc) etc, equation (4) acquires another formally tilding-symmetric matrix 
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The pair of the old variables a and a must be replaced by their tilding-invariant 
product = aa = ^ and its tilding-covariant complement p = a/a = 1/ p. In an 
opposite direction, whenever needed, we may re-construct a and a from the two 
quadratic relations a^ = p^ and a? — ^/ p, i.e., up to an inessential indeterminacy in 
sign. After we abbreviate 
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and postulate that Sq = 2 and cxi = 1, this procedure results in the conclusion that 
our recurrences may be re-written as the following sequence of the coupled pairs of 
the recurrent relations. 



/c = l,2,... . 



(13) 



One re-interprets eqs. (13) as the mere recurrent definition of the auxiliary sequence 
of functions of our auxiliary real variable ^, 



Ei=e-2, (72 = ^-3, E2 = f-4e + 2, .... 

We see that the functions Sfc(^) and ak+i{^) are, by their construction, both poly- 
nomials of the same degree k. 
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Our final change in the notation will prescribe replaced by ^ = 4a;^, with 
represented as = 2T2k{x) and with re-scaled into = T2k-i{x)/x. We notice 
that our recurrences become simpler in the new notation but what is more important 
is that after such a transformation, our new polynomials T^^x) coincide precisely with 
the classical orthogonal Chebyshev polynomials of the first kind [12], 

To{x) = 1, Ti{x) = X, T2{x) = 2x^-1, .... (14) 

In this sense, our N — 3 recurrences are solved exactly in closed form. 



3.3 The two tilding-covariant variables at N = 4 



Using the abbreviations a — si and a — Sq etc., let us interpret the N — A QE 
recurrences (4) as a tilding-symmetric problem 
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The first and the last line may be dropped as defining merely p = —t/a andp = —i/a. 
Prom the next two outer lines we may express t in terms of i or vice versa. Once 
we normalize t = 1 and rc-parametrize — P = P{b) = + 2ha and, in parallel, 
Sq_i = (3 = (3{h) = (5^ + 2 6(5, the vanishing of the related secular two-by-two 
determinant may be re-read as the constraint 66=1. The rest of eq. (15) reads 
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Its first row expresses 7 as a function of a and b and a and b, 

^yba/a — Q — P{b) ba + aa — 3 . 



(17) 



Fortunately, the quantity Q may be read as a polynomial in the mere two new, 
auxiliary variables Z = a a = Z and Y = ab = Z/Y. In the light of the preceding 
subsection, this function Q — Q{Z, Y) — ZY + 3Z — 3oi two variables will obviously 



7 



play the role of a generalized Hermite polynomial. We may only regret that it is not 
a tilding-invariant function anymore. 

The second row of eq. (16) must be multiplied by Fa to give the next-order 
polynomial in the same two variables, 

R^5ha^Y/a^ ZYQ{Z,Y) + Z^Y + 2Z'^ - ZY ^ R{Z, Y) . (18) 

This defines the new quantity (= a rescaled 5) and we may proceed to the third row 
multiplied by Y a'^, 

S = ebd? Y/a = Z R{Z, Y) + Z^ Q{Z,Y) - Z^Y - 2 Z^ = S{Z, Y) . (19) 

After the multiphcation by Y'^ the fourth row reads 

T = Cba^ Y^a = ZY S{Z, Y) + Z^ R{Z, Y)-YZ^ Q{Z, Y) (20) 

with the next factor Y"^ giving the next, fifth row 

U^rjba^ Y^/a = Z T{Z, Y) + Z^ S{Z, Y) - Z^ R{Z, Y) (21) 

etc. Step by step we construct, in this manner, the two sequences of functions denoted 
as Pn{Y, Z) and Qn{Y, Z) and defined by the pair of the coupled recurrences, 

Pn+l — YZ Qn + Z"^ Pji — Y Z^ Qn-l-i 

Qn+1 — Z Pn+l + Z^ Qn — Z^ P„, (22) 

n = 0, 1,... 

from the initial values Q-i = ^/Z, Pq — Y+2 and Qo — Q — Y Z+3 Z — 3 generating 
R^Pi etc. 

In a way paralleling the previous = 3 case, we might slightly modify the 
functions and define Pn+i = VY W2n+i while Qn+i = W^2n+2- It is easy to verify that 
we can now use just the single common recurrence 

Wn+l ^WZWn + Z^ Wn-1 - W Z^ Wn-2 , /,,oN 

n = 0,l,... ^^"^^ 

with the merely slightly modified initialization by W-2 — Q-i — 1/-^, W-i — 
Po/VY ^ VY + 2/Vy and Wo = Qo = Q ^ Y Z + 3Z - 3. We may see the 
clear parallels with the previous = 3 case, noticing that the polynomials W^n and 
Wsn-i are both divisible by Z^" while H^3n-2 is only divisible by We shall 

skip the further technical analyses of this sort here. 
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4 Matching and secular polynomials at N = 3 



For the sake of simplicity, let us only pay attention to the choice of = 3. Then, 
the knowledge of the closed form of the polynomials T,k{0 ^^<i crfe+i(0 enables us to 
define the explicit values of all our coupling constants as functions of the mere two 
parameters a and d entering ^ — aa and p — a/d, 

ai = a = ap°(7i(^), A = ^ = aa = 6 = ap(72(0, 

^ = 5 = p'S2(0, a3 = c = apV3(0, ••• • (24) 

One could also have constructed this general solution of our recurrences (10) in an 
opposite, upward direction. For this purpose, it suffices when all the above formulae 
are modified by a consequent application of the tilding operation. 

We have seen that the N — 3 case operates with two unknowns. At the same time, 
the set of recurrences (10) contains precisely two redundant items. In one extreme 
example we may read whole this set as a sequence of definitions of j3 = S2 = S2(si, Sg), 
7 = S3 = S3{si,Sq), . . . , = Sq+i{si,Sq) whcrc the last two lines are redundant 
since we already knew the outcome, viz., s^+i = 2 and Sq{si,Sq) = Sg. This may 
be understood as a source of our final pair of boundary conditions determining the 
QE-compatible values of the pair of the unknown parameters. 

In a way paralleling the previous N = 2 example, any other two lines of eq. (10) 
might be selected as boundary conditions. In contrast to the N = 2 example, almost 
all of the non-extreme choices of matching conditions would be preferable in practice, 
lowering the degree of the resulting secular polynomials. 

This observation deserves to be explained in more detail. Indeed, it makes sense to 
distinguish between the four possible selections of the optimal matching conditions. 

4.1 q = 4:K 

Whenever q = 4K where K = 1, 2, . . ., the above-mentioned recurrent construction 
may be started, simultaneously, at both the upper and lower end of eq. (10). Without 
any difficulties and using eq. (24), the recipe defines all the unknown quantities, i.e., 
the doublets of pairs of the values 

{aj, Aj), (dj, Aj), j^l,2,...,K. (25) 

The two middle lines of eq. (10) define the other two redundant functions (or "non- 
existent couplings") Qk+i and dx+i- This induces no real difficulty since the two 
parameters ai = a and di = d are not yet specified. The latter two definitions are 
not redundant, therefore, as they have to fix these initial values. 

The inspection of eq. (10) reveals that our symbol gk+i is an alternative name 
for another and well defined coupling A^- Similarly, the quantity dK+i is an "alias" 
for Ak- We determine the missing QE roots a and d via the two redundant equations 
Qk+i — Ak and dx+i — Ak or, in the notation of eq. (24), 

a gk+i ^T^kP^ , a ax+i p^ ^T,k p^ ■ 
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Their ratio reads p — 1 and gives the unique real root p = 1. Our first conclusion 
is that we must have a = a. The above two equations coincide and any of them 
represents our ultimate matching condition or constraint imposed upon ^ = a^, 

aaK+i{a^)^^K{a^), q ^ AK . (26) 

A sample of its roots may be found in Table 1 at g = 4 and q — 8. The inspection of 
the subsequent Table 2 reveals that with the further growth of K, the determination 
of these roots becomes purely numerical very quickly. 



4.2 q = 4:K + 2 

After a move to q — AK + 2 with K — (0, ) 1, 2, . . ., the previous recipe does not 
change too much. This time we define all the unknowns in a reversed order, 

{A,, a,+i) , (i„ a,+i) , J = 1, 2, . . . , . (27) 

Mutatis mutandis we find that the central part of eq. (10) defines the other "non- 
existent" couplings Ak+i and Ak+\ so that the doublet of equations Ak+i — ax+i 
and Ak+1 = ax+i leads to another set of the selfconsistency conditions, 

a gk+i = SftT+i p^^^ , a uk+i p^ = S^r+i p^^^ ■ 

Their ratio degenerates to the modified constraint = i with the same unique 

real root as above, a/a — p — 1. Both our innovated identities coincide, 

aaK+i{a^) = T.K+M), q^4K + 2 (28) 

and guarantee the desired matching. Their numerical aspects are sampled again in 
Table 1 (easily solvable cases at 7^ = 0, 1 and 2). The adjacent Table 2 complements 
this list and facilitates the determination of the explicit form of the secular equation 
(28) at all the integers K. 



4.3 q = 4K-\-l 

The subset of odd q = 4K + 1 with = (0, ) 1, 2, . . . requires a more careful analysis. 
Although we have the same complete list (27) of the definitions of the QE-fixed 
couplings as above, its last two items are defined twice, in two different ways. Their 
necessary compatibility represented by the relation a^+i = o,k+i or rather 

a ax+i p^ = a ax+i p^ 

implies that p^^"*"^ = 1 so that we must put a/a = p = \. In the light of this 
conclusion, the other two consequences A^+i = A^ and Aj^+i = A^ of the two 
other next-to-central rows of eq. (10) coincide and give the same ultimate matching 
rule 

J:k+i{0-^k{0, g = 4X + l. (29) 
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Its numerical performance appears illustrated by the corresponding subset of roots 
in Table 3. 

Marginally, let us note that for the specific exponents q = iK + 1, the secular 
polynomial may be re-written in the compact form 



(30) 





The secular polynomials J2m=o contain the {K + 1)— plets of coefficients 

/CW = (SP, dP) such that /C(o) = (1), /CW = (1, -1), /C(2) = (1, -3, 1), 

/C^^^ = (1, —5, 6, 1), etc. This rule parallels the even— g recipe of Table 2. 

4.4 q = 4:K + 3 

The last possible choice of the odd exponents q — AK + 3 (with K — (0, ) 1, 2, . . .) 

in the potentials V(g,n)(r) of cq. 2) leads to a routine completion of all the above 
analysis. A marginal modification of the list (25) is needed to specify all the necessary 
QE couplings, recurrently determined as functions of a and d only, 

{aj, Aj) , (dj, Aj), J = 1, 2, . . . , X + 1 . (31) 

Nevertheless, the results of the matching become slightly different this time. Al- 
though the first, central-line rule Ak+i = Ak+i prescribes merely 

^K+l P — P 

its consequence p^^^"^ — 1 admits the two alternative signs in the resulting a/d — 
p = ±1. Under this condition, the other two equations (in detail, aK-\-2 — S/f+i and 
its tilding-conjugate dK+2 — o^+i) coincide as well, giving the same final condition 

aK+2{i)=(JK+i{i), g = 4X + 3. (32) 

Curiously enough, this equation is the most easily solvable implicit definition of the 

QE roots ^ = a2 = 52 ^^^^ r^^^Ae 4). 

Even the shortest glimpse at the results of the factorization of the effective secular 
polynomial R^^''^\^) — ux+^iO —(^k+i{C) reveals that the sequence of the exponents 
q — AK + 3 might be viewed as the most privileged one. The search for its QE roots 
becomes by far the easiest. After we omit the roots ^ = and ^ = 4 as trivial, we 
encounter another unexpected and purely empirically observed symmetry. Indeed, 



the secular roots ai = y/^±n = y 2 ± s[f ' listed in Table 4 at the indices q = AK + 3 
appear to be of the very similar form. It is realy instructive to list a few sample 
distances S^f-l of al from their median = 2. We have = 0, = 1, Sifl = 
or ^/2, two values of = (-s/S ± l)/2, the three values of Slfl = 0, 1 and VS. One 
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may see that the full set of the secular roots oi — exhibits a weird regularity 
manifested by a reflection symmetry with respect to the center at = 2. All 
roots become tractable as certain quasi-conjugate pairs ^ = ^±n = 2 ± Stf-l. In this 
sense, the results hsted in Table 4 may be tentatively extrapolated to all the values 
of q. Indeed, once we omit the permanent pair of the minimal and maximal QE- 
compatiblity roots s = Si — ±2 as a trivial, we may use the auxiliary variable H^^l 
as defined by the relation Si = V2 ± at all the integers K. The inspection of 
Table 4 then reveals that all the complete sets of all the QE-compatiblity roots at 
all the hsted integers q — AK + 3 (i.e., at the set of X = 1, 2, 3, 4 and 5) coincide 
with the complete sets of roots of the Chebyshev polynomials of the second kind, 

^^(^^)"°' ^ = l'2,...,Xo (33) 

with the confirmed Kq = 5 at present. This means that in a way complementing 
the N > Nq extrapolations performed in our paper PI we may now tentatively 
extrapolate also the result (33) and conjecture that this equation determines all the 
closed and exact solutions of our N — 3 QE Magyari-Schrodinger eq. (4) also at all 
the larger integers q^AK + 3> AK^ + 3 = 23. 
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• Table 1. Real QE roots s = si at the first few even q for N = 3. 

• Table 2. Double Pascal triangle for coefficients in the reduced secular equations 
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• Table 3. Real QE roots s = Si at the first few odd q ior N — 3. 

• Table 4. Real QE roots s = Si at the first few q = 3 {mod 4) for = 3. 

Tables 



Table 1. Real QE roots s = si at the first few even q ioi N = 3. 



Q 


roots s — Si 


2 


2 




-1 






4 


2 


(V5- 


-l)/2 -(v/5 + l)/2 






6 


2 




(plus all three roots of + - 


- 2a 


-1) 


8 


2 




-1 (plus all three roots of 




-3a + 1) 


10 


2 


(plus 


all five roots of + — 4a^ — 


3a2 


+ 3a + 1) 
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Table 2. Double Pascal triangle for coefficients in the reduced secular equations 
J2k=o ^1 ^if — 0- This extrapolates Table 1 to all the even g's. 









coefficients 











1 


2 


3 


4 


5 ... 


Q 














2 


1 


1 










4 


— 1 


1 


1 








6 


-1 


-2 


1 


1 






8 


1 


-2 


-3 


1 


1 




10 


1 


3 


-3 


-4 


1 


1 


12 


-1 


3 


6 


-4 


-5 


1 ■-. 


14 


-1 


-4 


6 


10 


-5 


-6 ■-. 


16 


1 


-4 


-10 


10 


15 


-6 ■-. 



Table 3. Real QE roots s = Si at the first few odd q ios: N — 3. 



Q 






roots s 






1 
3 
5 
7 

9 
11 


2 
2 
2 
2 

2 
2 


-2 
-2 
1 

V2 


-1 

1 


-2 

-2 

-1 


^3+V5 2 

-v^ -2 



Table 4. Real QE roots s = si at the first few ^ = 3 {mod 4) for N — 3. 





roots s = Si 


3 


±2 






7 


±2 






11 


±2 


±V2±1 




15 


±2 




19 


±2 


±y2± (75 T l) /2 




23 


±2 


±^ 
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